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Abstract 

Nominal roll control laws were designed, 
implemented, and tested on an acroelastically-scalcd frec- 
to-roll wind-tunnel model of an advanced fighter 
configuration. The tests were performed in the NASA 
Langley Transonic Dynamics Tunnel. A parametric study 
of the nominal roll control system was conducted. This 
parametric study determined possible control system gain 
variations which yielded identical closed-loop stability 
(roll mode pole location) and identical roll response but 
different maximum control-surface deflections. 
Comparison of analytical predictions with wind-tunnel 
results was generally very good. 
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Abbreviations 


AFW Active Flexible Wing 

ARC Active Roll Control 

eg center of gravity 

cB decibel, 20 log(magnitudc) 

LEI leading-edge-inboard surface 

LEO leading-edge-outboard surface 
max maximum 

psf pounds per square foot 

TDT T ransonic Dynamics T unnel 

TEI trailing-edge-inboard surface 

TEO trailing-edge-oulboard surface 

A dot over a quantity indicates a derivative with respect to 
time. 


Consistent units have been used throughout the paper. 


Introduction 

The evolution of advanced fighters has required 
that the disciplines of aerodynamics, control systems, and 
structures be integrated into a unified aeroservoelastic 



technology which must be evaluated by sophisticated 
analytical methods and validated through the testing of 
wind-tunnel models. In support of this technology, 
Rockwell International Corporation has developed a 
concept it refers to as the Active Flexible Wing (AFW). 
This concept utilizes wing flexibility and multiple active 
control surfaces to vary the wing shape, resulting in 
improved performance and reduced weight. 

Under a joint Rockwell / Air Force / NASA 
program, using the AFW concept, Rockwell designed and 
built a scaled acroelastic wind-tunnel model which has 
been tested twice in the NASA Langley Transonic 
Dynamics Tunnel (TDT). A photograph of the wind- 
tunnel model mounted in the test section of the TDT is 
shown in figure 1. The AFW wind-tunnel model was 
mounted on a sting and attached to a large bearing 
arrangement in such a manner that it was free to roll. 
NASA synthesized Active Roll Control (ARC) control 
laws which were tested on the model during the second 
TDT entry. 

The purposes of this paper are to briefly describe 
the ARC control law design, to present and describe a 
NASA conducted control law parameterization study, and 
to present comparisons of predicted performance and 
predicted robustness with wind-tunnel test results. 


Wind-Tunnel Model 

The wind-tunnel model is an acroelastically- 
scalcd full-span model of an advanced fighter 
configuration. It has a fuselage and a low-aspcct-ratio 
wing with a span of approximately 9 feet. 

Slrucmrc 

The model consists of a "rigid" fuselage and a 
"flexible" wing. The fuselage contains aluminum 
stringers and bulkheads and is not scaled for flexibility. 
The wing box contains an aluminum honeycomb core and 
tailored plies of graphite-epoxy. The wing design 
permitted desired amounts of bending and twist as a 
function of aerodynamic load. The model was statically 
and dynamically scaled to represent a full-scale airplane 
with a wing span of roughly 50 feet. 

For acroelastic analysis purposes, the first 10 
calculated flexible antisymmetric mode shapes, 
frequencies, and generalized masses were provided to 
NASA by Rockwell International. Table I contains the 
natural frequencies and descriptions of these modes. All 
flexible modes were assumed to have a structural damping 
coefficient g of 0.03. 

Control Surfaces and Actuators 

The model has two leading-edge and two trailing- 
edge control surfaces on each wing panel. Each control 
surface has a chord of 25 percent of the local wing chord 


and a span of 28 percent of the wing semispan. Because 
the wind-tunnel model was aeroelastically scaled, the roll 
effectiveness of each pair of surfaces varied significantly 
with dynamic pressure. These effectivenesses were 
determined experimentally during the first TDT entry of 
this model. 

Each control surface is driven by a miniature 
rotary-vane electrohydraulic actuator which serves two 
functions: for constant inputs, it fixes the control -surface 
deflections relative to the wing; and for time-varying 
inputs, it provides control-surface motion in a manner 
dictated by the control law. Deflection limits are imposed 
on the various control surfaces to avoid exceeding actuator 
hinge moments and wing loads at the root. 

Instrumentation 

The model was instrumented with a force 
balance, accelerometers, strain-gage bridges, rotary 
variable differential transformers (RVDTs), a roll 
potentiometer, and a roll-rate gyro. Figure 2 contains a 
drawing of the model and illustrates the locations of the 
control surfaces and pertinent instrumentation. 


Wind Tunnel 

The NASA Langley Transonic Dynamics Tunnel 
(TDT) is a closed-circuit, continuous-flow tunnel which 
has a 16-ft. square test section with slots in all four walls. 
Mach number and dynamic pressure can be varied 
simultaneously, or independently, with either air or Freon 
as a test medium. All experimental tests of the present 
investigation were conducted in the TDT using a Freon 
test medium. 


Active Roll Control System 

At the time NASA was invited to participate in 
the AFW wind-tunnel tests the design objectives and the 
form (block diagram) of a roll control system, known as 
the Active Roll Control (ARC) System, had already been 
determined by Rockwell. 

Design Objectives 

The design objectives for the ARC system were 
as follows: 

Robustness 

0 + 6 dB gain margin 
0 +/- 45 degrees phase margin 
Performance 

0 Achieve 90 degrees roll angle in 
0.42 seconds 

The performance design objective is based on a scaled 
MIL-F-8785C (ref. 1), which specifies that the full-scale 
airplane be capable of rolling to 90 degrees in 1.1 seconds. 
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Block Diag ram 


The original form of the ARC system used a 
roll-rate gyro as the sensor. Rockwell's early analyses 
indicated that the ARC system would perform 
satisfactorily by using only the two most effective pairs 
of control surfaces at any test condition. Thus, it was 
necessary to include only two feedback loops and only two 
forward paths in the ARC system. 

In anticipation of a possible need to control 
flexible modes with the ARC system, NASA requested a 
modification to the original form of the block diagram, 
and a second-order filter was added to the then-existing 
system. The analog form of this second-order filter is 



s2 + 2 Cn a) n s+ “n 


A block diagram of the final form of the ARC 
system, including the NASA second-order filter, is shown 
in figure 3. The block diagram includes both digital and 
analog elements; the dashed line in the figure separates the 
digital and the analog portions of the closed-loop system. 

Control laws were implemented digitally with an 
Intel 80286 processor and 80287 co-processor. Tustin 
transformations (ref. 2) were used to generate the second- 
order filter in the feedback loop. The digital controller 
sampling rate was 200 Hertz. 

Design Conditions 

NASA designed nominal ARC control laws for 
two tunnel test conditions: a dynamic pressure of 150 psf 
at a Mach number of 0.9; and a dynamic pressure of 250 
psf also at a Mach number of 0.9. At the first condition 
the two most effective pairs of control surfaces were the 
trailing-edge-inboard pair and the trailing-edge-outboard 
pair; at the second condition, the trailing-edge-inboard pair 
and the leading-cdgc-outboard pair. 

Nominal Control Law Design 

For purposes of control law design, the 
acroclastic plant was represented by the rigid-body roll 
mode and the first 10 antisymmetric flexible modes of the 
acroclastic wind-tunnel model. Actuators were modeled 
analytically by zcroth-ordcr-over-third-order transfer 
functions. The coefficients of these transfer functions 
were obtained using parameter estimation techniques (ref. 
3) which matched both gain and phase from experimental 
transfer functions. 

Because the sampling rate of the digital controller 
was sufficiently high compared to the natural frequencies 
of key flexible modes, it was assumed that the digital 
controller was "approximately analog" and, therefore, the 
nominal control law was designed using classical analog 
techniques. The feedback gains Kj and K 2 , the gains 


K , and K _ in the forward path, and the constants C_ 
C 1 c 2 n 
and co n within the second-order filter were chosen by trial 

and error such that the performance and robustness design 
objectives were met analytically. 

Table II contains the results of the nominal ARC 
control-law design for the two test conditions. The left 
side of the table contains the gains and filter constants; the 
right side contains the predicted robustness and 
performance of the closed-loop system based on linear 
analysis (with no limits imposed on control-surface 
deflections). Nyquist plots (from which gain and phase 
margins were obtained) and time responses (from which 
time-to-roll-90-degrees were obtained) were computed 
using the ISAC and PADLOCS codes (refs. 4 and 5) with 
subsonic unsteady aerodynamics computed by the method 
of reference 6. The full equations of motion (rigid-body 
roll mode plus 10 flexible modes) and the full block 
diagram were present in the analysis. In addition, to 
approximate digital-to-analog conversion, a zero-order-hold 
element and a computational-delay effect were added to the 
analysis. The input used for the predicted performance 
calculations was a 0.04-second ramp-hold command with a 
magnitude of one radian per second. This input was 
chosen because the magnitude was sufficiently high to 
produce the desired predicted performance. Performance 
and robustness design objectives were met at both test 
conditions. 

Figure 4 contains two analytical Nyquist plots 
and illustrates (for an off-design condition) the potential 
effect of the second-order filter on the robustness of the 
closed-loop system. The analysis conditions for figure 4 
are a Mach number of 1 .15 and a dynamic pressure of 250 
psf. Again, the Nyquist plots were computed as above 
using the ISAC and PADLOCS codes but with supersonic 
unsteady aerodynamics computed by the method of 
reference 7. Figure 4(a) is for no second-order filter 
present and shows a clockwise encirclement of the minus- 
one point (at about 225 radians per second), corresponding 
to an instability of the fifth flexible mode. Figure 4(b) is 
for second-order filter present (with ^ = 0.8 and to n = 220 

radians per second), resulting in a very robust closed-loop 
system with over 10 dB gain margin and over 90 degrees 
phase margin. Comparison of the two Nyquist plots 
shows that stability was achieved by a 90-degree 
clockwise rotation and corresponding attenuation of the 
lobe associated with the previously unstable fifth flexible 
mode. 


Control Law Parameterization 
Background and Motivation 

The nominal ARC system has identical gains in 
each feedback loop and identical gains (but different from 
the values of the feedback gains) in each forward path. 
These equal values of gains result in certain closed-loop 
stability, certain closed-loop roll performance, and 
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(neglecting the effects of different actuator transfer 
functions) equal commanded control -surface deflections for 
the two pairs of surfaces involved. 

It was recognized (after making some simplifying 
assumptions) that there are an infinite number of 
combinations of feedback gains (Kj different from K 2 ) and 
an infinite number of combinations of forward-path gains 
( K . different from K_ ) which result in the same 
c | 2 

closed-loop stability and closed-loop roll performance but 
different commanded control-surface deflections for the 
surfaces involved. The consequence of this is that the 
deflections of one pair of surfaces may be "traded off 
against the deflections of the other pair with no change 
(loss or gain) of either stability or performance. 

With an eye toward actual airplane design, this 
trade off has beneficial implications in terms of wing 
loads, actuator sizing, and the ability of one pair of 
control surfaces to be used for multiple active control 
functions simultaneously. 

Development of Parametric Study 

The control law parameterization study was 
developed using the following guidelines. 

Feedback gains Kj and K 2 vary simultaneously 
(in general, Kj different from K 2 ) such that each is 

uniquely determined by a single parameter, K. Parameter K 
and gains Kj and K 2 arc continuously variable and any 

value of K (and, therefore, the corresponding unique pair of 
gains Kj and K 2 ) results in the same closed-loop stability 

as any other value of K. The relationship between 
parameter ic and gains K j and K 2 is described below in the 
subsection entitled "Constant Closed-Loop Stability." 

Forward-path gains K and K vary 
1 2 

simultaneously (in general, K different from K ) 

1 2 
such that each is uniquely determined by another single 

parameter, K_. Parameter tc„ and gains K. and K_ 
c c C 1 c 2 

are also continuously variable. For a given K and for a 

given roll-rate command, any value of k q (and, therefore, 

the corresponding unique pair of gains K and K ) 

C 1 C 2 

results in the same closed-loop roll performance as any 
other value of k £ . The relationship between parameter K c 

and gains K and K is described below in the 
1 2 

subsection entitled "Constant Closed-Loop Roll 
Performance." 


Choice of Stability and Performance 

For the purpose of this control law 
parameterization study the stated ARC robustness and 
performance design objectives cease to be important. It is 
required to show only that any desired robustness 
(stability) and any desired performance may be held 
constant for many values of K and for many values of k c . 

For convenience, the following stability and performance 
criteria were chosen for the parametric study: 

Stability 

0 Closed-loop roll-mode eigenvalue 
located at s = -26 on the negative 
real axis of the complex plane. 
Gain and phase margins not 
specified. 

Performance 

0 Achieve 90 degrees roll angle in 
1.5 seconds. 

The chosen stability represents a desired location of the 
roll-mode eigenvalue and has been scaled from the desired 
location for a full-scale airplane. The roll-rate command 
chosen to achieve this more conservative performance is a 
0.3-second ramp-hold command of the necessary 
magnitude. This command was chosen because it is a 
gradual, low-amplitude command which neither saturates 
the control surfaces nor excites flexible modes. 

Example Condition 

Whenever numerical examples arc required to 
illustrate the control law parameterization, the following 
test condition, referred to as the "example condition," is 
used: Mach number of 0.9; dynamic pressure of 250 psf. 
At the example condition the two most effective pairs of 
surfaces are the TEI and the LEO. To keep the 
illustrations general, the TEI surface is identified as 
Surface 1; the LEO surface as Surface 2. 

Problem Simplification 

The following approximations were made in 
order to simplify the mathematics while at the same time 
retaining the essential elements of control law 
parameterization: 

(1) From the ARC block diagram in figure 3, 
only the gains were retained. The following 
elements were neglected: stick-shaping 
filter, anti-aliasing filter, second-order filler, 
and actuator transfer functions. 

(2) From the aeroclastic plant, only the rigid- 
body roll mode was retained. All flexible 
modes were neglected. 

(3) The center of gravity of the wind-tunnel 
model was assumed to lie on the roll axis. 

Approximation (1) assumes that the neglected 
transfer functions are closely approximated by unity gain 
and zero phase over the frequency range of interest. The 
resulting simplified block diagram is shown in figure 5. 
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Approximation (2) deals with the aeroelastic 
plant (open-loop wind-tunnel model) only and assumes 
that the dynamics of the aeroelastic plant are closely 
approximated by the dynamics of the rigid-body roll mode. 
This approximation is valid if the frequency content of the 
commanded input is well below the frequency of the first 
flexible mode such that the flexible modes are not excited. 

The rigid-body roll equation for the open-loop 
wind-tunnel model is based on equation 4.15.8(b) in 
reference 8 and is given by 


I ♦ — TJSbt C $ +tnglsin$ =rJSbC 5 +l£bC 5 


( 1 ) 


P(s) 

P C ( S > 


A A 

K c L + K L 
1 S 1 2 5 2 
A A A 

s + L +K L + K L 
P 1 5 , 2 S 2 


W) 


A A A 

The quantities L„, L c , and L c arc 

P S 1 5 2 

constant at a given wind-tunnel test condition. Closed- 
loop stability may be deduced from the denominator of 
equation (4). For constant closed-loop stability (pole 
location) the quantity 


A A A 

L P +K 1 L 5i + K 2 L 6 2 


Equation (1) differs in two ways from the 
equation in reference 8 : control-surface rate derivatives 
have been neglected and the quantity mgl sin <}> has been 
added. This quantity is referred to as the "cg-offsct term” 
and is present because the center of gravity of the wind- 
tunnel model is a distance 1 below the roll axis of the 
model. The cg-offsct term causes equation (1) to be 
nonlinear^ The equation may be linearized by utilizing the 
small-angle approximation (substituting for sin <J>), 
resulting in the following 

1 xx* - ^ bt * S * + mgl 4>= qSbC , + cf.bC 1 6 2 < 2 > 


Approximation (3) assumes that the roll axis of 
the wind-tunnel model passes through its center of gravity 
(i.c., that the distance 1 is zero) as it would on an actual 
airplane. This approximation results in the following 
rigid-body roll equation 

I xx 4> -qSbt*C, 0 =qSbC, 8 , + qSbC , 5 2 (3) 


This assumption permits the derivation of 
control-law parameterization to proceed in a manner 
consistent with airplane equations of motion rather than in 
a manner consistent with wind-tunnel equations of 
motion. In a subsequent section of this paper showing 
comparisons of analytical predictions and experimental 
results the cg-offsct term will be added back into the 
analysis, as in equation ( 2 ). 

Constant Closed-Loop Stability 

Referring to the simplified block diagram in 
figure 5 and to equation (3) and recognizing that p = <j>, 
the following closed-loop transfer function is obtained 


must be constant, with value, say, s' (where, from an 
earlier section, s' = 26). The closed-loop eigenvalue will, 
therefore, be located at s = -s' on the negative real axis in 
the complex plane. By inspection, the following linear 
equation may be written for K 2 as a function of K j and 
results in this constant closed-loop stability 



All combinations of Kj and K 2 which satisfy equation (5) 
result in a closed-loop system whose eigenvalue is located 
at s = -s'. 


Figure 6 contains a plot of K 2 as a function of 
Kj and illustrates the correspondence between these gains 


and parameter k. The plot is a straight line lying in the 
second, third, and fourth quadrants of the Kj-K 2 plane. 

A 

L8 , 

The line has a slope of - 7 ; — , a K. -intercept of 

1 * 



, and a K 2 -interccpt of 



For the 


portion of the line in the third quadrant, parameter K varies 
linearly, as shown in the figure, with values between zero 
and one. The following expressions define gains Kj and 

K 2 as functions of parameter tc 


K 


1 



K 


( 6 ) 
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s’ -L 

^ 2 = a (1 " K) (7) 

\ 

When K is equal to zero, gain Kj is zero (Loop 1 
is open) and gain K 2 alone is "holding" the closed-loop 

eigenvalue at s = -s'. Values of k less than zero (values of 
K j greater than zero) corresponds to a stable closed-loop 
system in which Loop 1, by itself, is unstable. 

When K is equal to one, gain K 2 is zero (Loop 2 
is open) and gain Kj alone is "holding" the closed-loop 

eigenvalue at s = -s'. Values of tc greater than one (values 
of K 2 greater than zero) corresponds to a stable closed- 
loop system in which Loop 2, by itself, is unstable. 

The dashed line in figure 6 represents K 2 = Kj. 
The point of intersection of the dashed line with the solid 
line yields the value of k (= 0.76 for the example 
condition) which corresponds to equal values of feedback 
gain. 


Figure 7 contains on the left, the closed-loop 
eigenvalue in the complex plane and on the right, a 
Nyquist plot with the loop broken at the plant output. 
For these plots K = 0.76, but control law parameterization 
guarantees that for any k the closed-loop eigenvalue and 
the Nyquist plot will be identical to those in figure 7. 

Cgflgtao lQ psctl-LQQp RpH Performance 

For present purposes, "roll performance" is 
understood to be roll angle and roll rate as functions of 
time. Assuming a given roll-rate command and assuming 
having fixed the denominator of equation (4) by the 
method of the previous section, constant closed-loop roll 
performance may be obtained by fixing the numerator of 
equation (4) in a similar manner. If the value of the 
numerator is arbitrarily set to the value which results 
when K and K arc both equal to K ', the 
C 1 2 

following linear equation, for K as a function of 

C 2 

K , may be written 



All combinations of K and K which satisfy 
Cl C 2 

equation ( 8 ) result in a closed-loop system whose closed- 
loop roll performance is constant. 


Figure 8 contains a plot of K as a function 

c 2 

of K c and illustrates the correspondence between these 


gains and parameter < c . The plot is a straight line lying 
in the first, second, and fourth quadrants of the K c - 


A 

L, 


1 


K plane. The line has a slope of - 7 — , a K 


intercept of 


/ A A 

K ’c L S +L * 

l 8 1 8 2 

A 

L_ 


, and a K c -intercept of 



For the portion of the line in the 


first quadrant parameter k c varies linearly, as shown, with 
values between zero and one. The following expressions 
define gains K c and K c as functions of parameter k c 



When K c is equal to zero, gain K is zero; 

C 1 

when k c is equal to one, gain K c is zero. The dashed 

2 

line represents K = K , and the point of 
2 c l 

intersection of the dashed line with the solid line yields 
the value of K c (= 0.76 for the example condition) which 
corresponds to equal values of gain in the forward path. 

Figure 9 contains analytical plots of the 0.3- 
sccond ramp-hold command and the resulting roll-angle 
and roll -rate responses for k = 0.76 and k c = 0.76. The 
magnitude of the command was chosen such that a roll 
angle of 90 degrees is achieved at time 1 .5 seconds. These 
time histories were obtained by applying the command to 
the closed-loop system represented by the simplified block 
diagram in figure 5. For these plots k = 0.76 and k £ = 
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0.76, but control law parameterization guarantees that for 
any k and for any K_ the roll-angle and roll-rate responses 
will be identical to those in figure 9. 

Control-Surface Deflections as Functions of K and 

Control-surface-deflection time histories vary as 
functions of parameters K and k £ . Figure 10 contains time 
histories of Surface 2 deflecting in response to the 0.3- 
sccond roll-rate command for three values of K (at K c = 
0.76). The maximum (absolute) value of each of these 
time histories is indicated on the plots with open circles. 
Figure 1 1 contains time histories of Surface 2 deflecting 
in response to the 0.3-second roll-rate command for three 
values of k c (at K = 0.76). The maximum (absolute) 
value of each of these time histories is indicated with 
closed circles. It can be seen from figures 10 and 1 1 that 
the character and magnitude of the deflection time histories 
can change significantly as k and k £ vary. 

It can be cumbersome and confusing to try to 
understand the variation of entire time histories with 
respect to even one parameter. A better approach would 
be to try to understand the variation of a single quantity 
related to each time history. Such an approach is taken 
here. The measure of how these time histories vary is the 
absolute value of the maximum deflection obtained in 
response to the 0.3-second roll-rate command. This 
section presents the variation of these maximum 
deflections for values of k between zero and one and for 
values of * c between zero and one. 

Figure 12(a) contains an analytically-predicted 
contour plot of the maximum values of Surface 1 

deflection (6 ) as functions of k and k„. It was 

1 w 
max 

obtained by sampling many time histories of 8 j as K and 
K c were varied. Parameters K and k £ were each varied from 
0.00 to 1.00 in increments of 0.05, for a total of 441 
combinations of k and k c . The plot is a "valley" with the 
locus of minima running diagonally. The contour interval 
is 1 degree. The minimum point on the plot is zero and 
occurs in the lower left comer (at k = k £ = 0), which 
corresponds to K i = K =0. From the simplified block 
C 1 

diagram in figure 5, K j = K c =0 means that Surface 1 

receives no signals to deflect: none from the roll-rate 
command; none from the roll -rate feedback. The 
deflection of Surface 1 is therefore zero and Surface 2 does 
all the work. 

Figure 12(b) contains an analytically-predicted 
contour plot of the maximum values of Surface 2 
deflection ( 5 )as functions of K and K c . The plot is 
2 max 


also a "valley" with the locus of minima running 
diagonally. Again, the contour interval is 1 degree. The 
open and closed circles correspond to those from figures 
10 and 1 1. The minimum point on the plot is zero and 
occurs in the upper right comer (at k = K c = 1), which 

corresponds to K 0 = K =0. From the simplified block 
L C 2 

diagram in figure 5, K-, = K =0 means that Surface 2 

C 2 

receives no signals to deflect: none from the roll-rate 
command; none from the roll-rate feedback. The 
deflection of Surface 2 is therefore zero and Surface 1 does 
all the work. 

Trading Off Control-Surface Deflections 
Against Each Other 

By proper choice of the parameters K and the 

maximum deflection of one control surface may be traded 
off against the maximum deflection of the other. As an 
illustration, figure 13 contains plots of 8 and 

* max 

8 as functions of k, for tc = 0.76. The curves in 
2 c 
max 

figure 13 correspond to taking "slices" (at this value of K c ) 
through the contour plots in figures 12(a) and 12(b). In 
figure 13, deflection 8 reaches a minimum value at 
2 max 

K = 0.62, which corresponds to a point at the bottom of 
the "valley" of Figure 12(b). In figure 13, however, 
deflection 8 decreases over the entire k range 
1 max 

without reaching a minimum. (It can be seen from figure 
1 2(a) that, at this value of k c , the bottom of the "valley" 

occurs at a value of k greater than one.) From Figure 1 3 it 

can been seen that 8 and 8 are equal at two 

12 
max max 

values of k. In the range of k between the two values, 

8 is less than 8 ; outside this range of k, 

max 1 max 

8 is greater than 8 . At k = 0.76, not only 

L max 1 max 

are the maximum values equal, but the entire time 
histories are identical. 

Thus, if it were necessary (because of, say, the 
resulting reductions of wing loads) to minimize the 
maximum deflection of Surface 2 (while at the same time 
maintaining stability and performance constant), the 
control law parameterization provides the value of K (and, 
therefore, the gain combination, Kj and K^) necessary to 
accomplish this task. 
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Wind-Tunnel Results and Comparison with Analysis 
Frequency-Response Te st Technique 

Nyquist frequency-response plots were obtained 
experimentally with the wings level and with the wind- 
tunnel model free to roll. In the ARC system the 
feedback path was broken at the plant output and an 
analog sinusoidal signal (whose frequency varied 
logarithmically between 2 and 20 Hertz over a span of 
about two minutes) was inserted at the break. This signal 
and the analog output signal from the roll-rate gyro were 
recorded and later processed by an HP-5420 signal 
analyzer, yielding the experimental Nyquist plots. 

Deflection limits were imposed on the control 
surfaces at all times. Because of this imposition, care was 
taken to insure that the magnitude of the input sinusoidal 
signal was small enough so that the control surfaces 
would never saturate at any time during the Nyquist tests. 

Time-Response Test Technique 

Commanded roll maneuvers were performed 
using the digital control computer. The wind-tunnel 
model was initially held in place with the left wing 45 
degrees down by a roll-trim system also residing in the 
digital computer. At the initiation of the maneuver the 
computer disengaged the roll-trim system and engaged the 
roll-control system. Referring to the block diagram in 
figure 3, the computer then generated the 0.3-second ramp- 
hold roll-rate command signal. This signal passed 
through the stick-shaping filter and forward-path gains 
through the digital-to-analog converters and into the 
actuators. In response to the deflecting control surfaces, 
the model rolled through the wings-levcl position to a 
position with the left wing 45 degrees up, for a total 
incremental roll angle of 90 degrees. At this point the 
computer terminated the maneuver by re-engaging the roll- 
trim system to hold the model in the left-wing-up 
position. 

The 45-dcgrcc start and stop positions were 
chosen because of the small angle approximation made in 
the cg-offsct term of the rigid-body roll equation (eq. 2). 
The difference between <!> and sin (j> is only about 10 
percent when <J» is 45 degrees, and this difference 
.approaches zero as <J> approaches zero. Thus, this choice 
of start and stop positions minimizes the error incurred 
because of the approximation. 

Correction to Experimental Data 

About midway through the three-week wind- 
tunnel test, it was determined from comparison of 
integrated roll rate with roll-angle measurement that the 
roll rate gyro output was 25 percent too high. At that 
lime, and for all subsequent runs, a factor of 0.8 was 
applied to the output of the roll-rate gyro to account for 
the discrepancy. All experimental data at the first test 
condition (M = 0.9, q = 150 psf) was obtained before the 


0.8 factor was applied to the output of the gyro. All 
experimental data at the second test condition (M = 0.9, q 
= 250 psf) was obtained after the 0.8 factor was applied to 
the output of the gyro. Therefore, for consistency in 
comparing wind-tunnel data with analysis, any 
comparisons at the first test condition have been rescaled 
to include the 0.8 factor. 

Nominal Control Law 

Figure 14 contains a comparison of analytical 
and experimental Nyquist plots for the nominal ARC 
control law at a Mach number of 0.9 and a dynamic 
pressure of 150 psf. The analytical Nyquist plot was 
obtained using, again, the ISAC and PADLOCS codes 
with the full equations of motion (rigid-body roll mode 
plus ten flexible modes) and with the full block diagram 
(all filters as well as the zero-order-hold element and the 
computational-delay effect) present in the analysis. The 
cg-offset term was not present in the analysis. The 
experimental Nyquist plot was obtained in the manner 
described above. It can be seen from figure 14 that the 
analysis correctly predicted the general shape of the plot as 
well as the relationships between the first two flexible 
modes. The analysis predicted the phase margin well but 
was unconservative (by about 3 dB) in its prediction of the 
gain margin. 

Table III contains a summary of the analytical 
and experimental gain and phase margins for both test 
conditions. The gain-margin design objective was met at 
one test condition; the phase-margin design objective was 
met at both. 

Control Law Parameterization 

From an earlier section of this paper the control 
law parameterization was constructed to yield identical 
Nyquist plots for all values of k. Because the derivation 
of control law parameterization was based on a single- 
degree-of-freedom system this requirement holds for a 
single-degree-of -freedom system. However, for a many- 
degree-of-freedom system (such as the AFW wind-tunnel 
model), which is closely approximated (at low frequencies) 
by a single-degree-of-freedom system, the requirement 
should hold only approximately. 

To investigate the validity of the 
parameterization Nyquist plots were obtained analytically 
and experimentally for several values of k between zero 
and one, and gain margins and phase margins were taken 
from the plots. The analytical Nyquist plots were, again, 
obtained with the full equations of motion and with the 
full block diagram present in the analysis. 

Figure 15 contains a comparison of analytical 
and experimental gain and phase margins as a function of 
parameter k at a Mach number of 0.9 and a dynamic 
pressure of 150 psf. It can be seen from the figure that 
analytical gain and phase margins and experimental gain 
and phase margins are approximately constant as a 
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function of parameter tc. Gain margins only vary by 
about +/- 1 dB; phase margins vary by about +/- 5 
degrees, (Although not shown, the Nyquist plots are also 
approximately constant as a function of parameter k.) 
This trend verifies the control law parameterization. 

To investigate the trading off of control surface 
deflections as predicted by the control law 
parameterization, time responses were obtained 
analytically and experimentally for several values of tc and 
kt c between zero and one. The analytical time responses 

were obtained with the rigid-body roll-mode 
approximation to the full equations of motion (cg-offset 
term present — cq. (2)) and with the simplified block 
diagram (fig. 6). The roll-rate command, analytically and 
experimentally, was the 0.3-second ramp-hold signal with 
the same magnitude. 

Figure 16(a) contains a comparison of maximum 
(absolute) values of deflections 8j (TEI) and §2 (LEO) as 
a function of k (for k. = 0.76) at a Mach number of 0.9 
and a dynamic pressure of 250 psf; figure 16(b) contains a 
similar comparison as a function of < c (for K = 0.76) . 
The solid and dashed lines correspond to the analytical 
predictions of 8 and 8 , respectively; the open 

* max ^ max 

and closed symbols correspond to the experimental results. 
From figure 16 it can be seen that the analysis correctly 
predicts the behavior of 8 and 8 as functions 
1 max ^max 

of K and k c , and, therefore, verifies the trade-off predicted 
by the control law parameterization. 

Concluding Remarks 

This paper (1) outlines the design and 
implementation of nominal control laws for active roll 
control, (2) describes a control law parameterization study, 
and (3) presents a comparison of wind-tunnel results and 
analytical predictions for an aeroelastically-scaled wind- 
tunnel model of an advanced fighter configuration. The 
model was free to roll and was tested in the NASA 
Langley Transonic Dynamics Tunnel. The nominal 
control laws met the gain-margin design objective (+ 6 
dB) at one test condition and met the phase-margin design 
objective (+/- 45 degrees) at both test conditions. 
Comparison of analytical predictions with wind-tunnel 
results was generally very good and verified both the 
robustness and the performance requirements for the 
control law parameterization study. 
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Table I. - Structural Modes of Wind-Tunnel Model 



Description 

Frequency, Hz 

1 

Sting 1 st bending 

7.18 


Wing 1 st bending 

12.83 

H 

Fuselage yaw 
Wing 2nd bending 

16.60 

34.17 


Wing 1 st torsion 

35.05 


Wing/fuselage 

38.59 

IBsP 

Wing 3rd bending 

47.96 


Wing bending/torsion 

51.10 

9 

Wing 2nd torsion 

53.03 

10 

Wing/fuselage 

56.99 


Figure 1. - AFW wind-tunnel model in NASA Langley 
Transonic Dynamics Tunnel. 



Table II. - Nominal Control Laws for Active Roll Control 
System Mach = 0.9 


q 


k 2 


“n 

Gain margin 

Phase margin 

Time 
to 90 

(psf) 

K 1 


dB 

@o> 

Deg. 

@co 

(fiS) 

(sec) 

150 

0.2 

(TEI) 

0.2 

(TEO) 

0.8 

220 

8.3 

60 

72 

21 

0.28 

250 

0.2 

(TEI) 

0.2 

(LEO) 

0.8 

220 

9.4 

63 

76 

23 

0.30 


6 DOF force balance 



1 Roll set for 
roll angle 

11 Accelerometers 
5 per wing 
1 on fuselage 


8 - RVTD for surface 
position 

1 per control surface 


16 Strain gages 

1 Roll rate gyro • 

Figure 2. - Wind-tunnel model instrumentation. 


Table III. - Comparison of Robustness Characteristics for 
Nominal Control Laws 



Gain margins 

Phase margins 

Dynamic 

Analysis 

Experiment 

Analysis 

Experiment 

pressure, 

psf 

dB 

ii 

dB 

4 1 

Deg. 

life 

Deg. 

ii 

150 

8.3 

60 

5.4 

57 

72 

21 

68 

25 

250 

9.4 

63 

6.5 

55 

76 

23 

63 

25 



Figure 3. - Block diagram of active roll control system. 


10 


ORIGINAL PAGE IS 
OE POOR QUALITY 



















(a) Without second-order filter 
(unstable) 


(b) With second-order filter 
(stable) 


Figure 4. - Effect of second-order filter on Nyquist plots. Mach =1.15; dynamic pressure = 250 psf. 



Figure 5. - Simplified block diagram of active roll 
control system. 


Figure 6. - Feedback gains which result in constant 
closed-loop stability. 
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(a) Eigenvalue (b) Nyquistplot 

Figure 7. - Closed-loop eigenvalue and associated Nyquist plot for constant closed-loopstabilily. 



* 


4-c/V 

Figure 8. - Forward-path gains which result in constant closed-loop roll performance. 
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(a) Roll-rate command 




0 


1 

Time, sec 

(c) Roll-angle response 



Time, sec 
(a) k = 0.56 



Time, sec 
(b) k = 0.76 



Time, sec 
(c) k = 0.96 


Figure 9. - Ramp-hold command and constant-loop roll 
performance. 


Figure 10. - Time histories of Surface 2 (LEO) deflections as a 


function of K. k = 0.76. 

C 
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Figure 13. - "Slices" through contour plots for constant k . 



Figure 14. - Comparison of analytical and experimental Nyquist plots. Nominal control law; 
Mach = 0.9; dynamic pressure = 150 psf. 
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10 

8 

Gain 6 

margin, 
dS 4 

2 

0 

90 


Phase Rn 
margin, 
degrees 

30 


0 .5 .69 .9 1.0 

K 

Figure 15. - Comparison of analytical and experimental gain and phase margins as functions of K. Mach = 0.9; 

dynamic pressure =150 psf. 



0 

0.5 

1.0 

0 0.5 

1.0 


K 





(a) k variation 


(b) k c variation 



*c = 0.76 


k = 0.76 



Figure 16. - Comparison of analytical and experimental parameterization results. Mach = 0.9; 

dynamic pressure = 250 psf. 
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